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We compute the photon creation inside a perfectly conducting, three dimensional oscillating cav- 
ity, taking the polarization of the electromagnetic field into account. As the boundary conditions 
for this field are both of Dirichlet and (generalized) Neumann type, we analyze as a preliminary 
step the dynamical Casimir effect for a scalar field satisfying generalized Neumann boundary condi- 
tions. We show that particle production is enhanced with respect to the case of Dirichlet boundary 
conditions. Then we consider the transverse electric and transverse magnetic polarizations of the 
electromagnetic field. For resonant frequencies, the total number of photons grows exponentially in 
time for both polarizations, the rate being greater for transverse magnetic modes. 



I. INTRODUCTION 



' The existence of an attractive force between two uncharged, perfectly conducting parallel plates was predicted by 
T^J- , Casimir in 1948 jl| and has recently been measured at the 15% precision level using state-of-the-art cantilevers [|J. 

A similar force between a conducting plane and a sphere has also been measured with progresively higher precision 
in the last years using torsion balances ||, atomic force microscopes [Q, and capacitance bridges f|||, with the 
latter reference showing the relevance of Casimir forces in nanotechnology. For a recent review of experimental and 
theoretical developments see 0. 

The dynamical effect consists in the generation of photons due to the instability of the vacuum state of the electro- 
magnetic field in the presence of time-dependent boundaries. In the literature it is referred to as dynamical Casimir 
effect ||] or motion- induced radiation ||. The dynamical effect has been recently reviewed in Ref. p0| . Up to now 
no concrete experiment has been performed to confirm this photon generation, but an experimental verification is 
not out of reach. From the theoretical point of view it is widely accepted that the most favorable configuration in 
order to observe the phenomenon is a vibrating cavity in which it is possible to produce resonant effects between the 
^ \ mechanical and field oscillations. 

Many previous papers have focused their attention in the scalar field quantization within a one dimensional cavity 
|Tl],[l2| . Recently, we analyzed in detail the case of a three dimensional cavity (ll| , but still considering a scalar field 
(in this and other ]l4| previous works Dirichlet boundary conditions are assumed). The main difference between one 
and three dimensional cavities is that, while in one dimension the cavity's frequency spectrum is equidistant and leads 
to strong intermode interactions, in three dimensions the spectrum is in general nonequidistant, and only a few modes 
may be coupled. The relevance of this coupling has been pointed out only recently (see fU^Ts! ) . The aim of this 
paper is to extend the results of Ref. |l3| to the realistic case of the electromagnetic field, properly taking into account 
the polarization of the different modes (transverse electric (TE) and transverse magnetic (TM) polarizations). 

As we will see in Sec. II, the electromagnetic field involves both Dirichlet and (generalized) Neumann boundary 
conditions. For this reason, it is of interest to analyze the case of a massless scalar field satisfying this latter type of 
boundary conditions, which we do in Sec. III. Assuming a resonant vibration of the cavity and using multiple scale 
analysis (MSA) |H| we will show that the number of particles produced is much larger than for Dirichlet boundary 
conditions. We study in detail the resonant case in which the cavity oscillates at twice the frequency of some field 
mode. In Sec. IV we show that TE modes of the electromagnetic field behave as a scalar field with Dirichlet boundary 
conditions, while TM modes are analogous to the scalar case with Neumann boundary conditions. Section V contains 
our main conclusions. 
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II. THE BOUNDARY CONDITIONS 



We consider a rectangular cavity formed by perfectly reflecting walls with dimensions L x ,L y , and L z . The wall 
placed at x — L x is at rest for t < and begins to move following a given trajectory, L x (t), at t = 0. We assume this 
trajectory as prescribed for the problem (not a dynamical variable) and that it works as a time-dependent boundary 
condition for the field. Moreover, we will assume a non relativistic motion of the wall with L x (t) = Lq(1 + e/(i)) with 
e<1 and f(t) a bounded function. We use units K = c = 1. 

Let us consider the problem of finding the electromagnetic field inside the cavity in terms of the four vector potential 
= (tp, A). In the Coulomb gauge V • A = the scalar potential tp vanishes and the vector potential satisfies the 
wave equation DA = 0. For the static walls, the boundary conditions are the usual ones 

E||=0 ; B ± =0, (1) 

where || and _L respectively denote the components of the field parallel and perpendicular to the wall. 

On the moving wall, these boundary conditions must be imposed in a Lorentz frame in which the mirror is instanta- 
neously at rest. As the mirror moves in the x direction, it will be convenient to decompose the electromagnetic fields 
into TE and TM modes with respect to the x axis. The TE fields are defined as the solutions to Maxwell equations 
with E( TE ) • x = 0. Analogously, the TM fields satisfy B^™) • x = 0. A general solution for the electromagnetic field 
inside the cavity is a superposition of TE and TM fields. 

Hi* 



It is useful to introduce a different vector potential for each polarization through the equations 1 17 



E< TE ) = -cW TE) ; B( TE ) = V x A( TE >; 
B(™)=d t A(™) ; e(™) = VxA(™1 



(2) 
(3) 



Both potentials satisfy the Coulomb gauge and have vanishing x component. As A • x = and tp = 0, the vector 
potentials are invariant under a boost in the x-direction. The same is true for the Coulomb gauge. In terms of 
these potentials, the boundary conditions are relatively simple Let us denote by S the laboratory frame and 

by S' the instantaneous comoving frame. In 5" the TE vector potential satisfies Dirichlet boundary conditions 
A'( TE )(V = 0,y', z',t') = 0. Therefore, on the moving mirror, 

A^ E \x = L x (t),y,z,t)=0. (4) 

On the other hand, the TM vector potential satisfies 

n^'d^A'^ix' = 0,y',z',t')=0, (5) 

where n M = (0, 1, 0, 0). As a consequence, for a non relativistic motion of the mirror, 

n^.Af™)^ = L x (t),y )Z) t) = (d x + L x (t)d t )A™(x = L x (t),y,z,t) = 0, (6) 

i.e., a "generalized" Neumann boundary condition with n M = [L x , 1,0,0). On the static mirrors the boundary 
conditions for the TE vector potential are given by 

A( TE )(x = 0,y,z,t)-O, 

4 TE > (x, y, {z = 0, L 2 }, t) = 4 TE) {x, {y = 0, L y }, z,t) = Q. (7) 
For the TM potential we have 

9*4™) (a: = 0, y, z, t) = d x A™ ( x = 0) y )Z)t ) = , 

4 TM) (*. {y = °> L y}' z ' *) = d v A ? M) (*> {y = o,L y },z,t) = o, 

4™> (x, y, {z = 0,L z },t) = 9,4™) (x, y, {z = 0, L 2 }, t) - 0. (8) 

From these boundary conditions it is clear that the behavior of each component of the TE vector field is related to the 
problem of a scalar field subjected to Dirichlet boundary conditions. For the TM vector field it is necessary to deal 
with the generalized Neumann boundary conditions given in Eq.(|^). The former problem was extensively studied in 
our previous paper p3|, and the latter one will be treated in the following section. 
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III. SCALAR FIELD WITH NEUMANN BOUNDARY CONDITIONS 



Let us consider the problem of a massless scalar field </>(x, i) satisfying the wave equation □</> = and (generalized) 
Neumann boundary conditions on all surfaces of the cavity. In the comoving frame the Neumann boundary condition 
is d^4>' = 0. In the laboratory frame, this condition becomes n^d^cj) = 0, where n M = (L x , 1,0,0). Therefore we 
have 

d x (j){x = 0,2/, z,t) = ; (d x + L x d t )(f>(x = L x (t),y,z,t) = 0; 
d y (t>{x, {y = 0, L y }, z,t)=0; d^x, y, {z = 0, L z }, t) = 0. (9) 



A. Instantaneous basis 



The Fourier expansion of the field for an arbitrary moment of time can be written in terms of creation and 
annihilation operators as 

0(x,t) =^a>„(x,t)+H.c, (10) 

n 

where the mode functions it n (x, t) form a complete orthonormal set of solutions of the wave equation with Neumann 
boundary conditions. When t < (static cavity) each field mode is determined by three non-negative integers n x ,n y 
and n z . Namely 



u n (x, t < 0) = 





with the shorthand n — (n x , n v ,n z ). In order to satisfy the boundary conditions for t > it is useful to expand the 
mode functions in Eq. (Un) with respect to an instantaneous basis. If the boundary condition on the moving mirror 
were the instantaneous Neumann condition d x (fi(x = L x (t),y, z,t) = 0, the trivial choice for the instantaneous basis 
would be 




However, as the generalized Neumann condition in Eq.(g) involves the time derivative of the field, the situation is 
more complex. 

Let us consider new variables (77, £) in the (t,x) plane. We define the line n = const as a slight modification of 
the corresponding t = const line, in such a way that it is orthogonal to the worldlines of the mirrors at x = and 
x = L x (t) (see Figure). The variable £ is defined as the distance, on the line r\ = const, from x — to x. In these 
coordinates, the generalized boundary condition on the two mirrors becomes d^(f>(^,y, z,rf) = both at £ = and 
at £ = l(r)), where l(rj) is the value of the coordinate £ on the moving mirror. Therefore, an instantaneous basis to 
describe the field is 
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FIG. 1. Worldlines of the mirrors along the x direction. The line rj = const is orthogonal to the worldlines at x = and 
x = L x (t), and the coordinate £ measures the distance from the static mirror along the r\ — const line. 



To find a concrete form for the new coordinates we write rj = t + g{x, i). Therefore £ is given by 

i / dx'Jl 



g' 2 (x',t) 



(l + g{x>,t)? 



(15) 



At this point it is important to note that, since we are considering motions of the wall that are small (0(e)) 
deviations from the initial static position, terms of order 0(e 2 ) or higher will be neglected in what follows. Moreover, 
g(x, t) = 0(e), £ = x + 0(e 2 ), and l(rf) = L x (t) + 0(e 2 ). With this in mind, it is easy to show that, in order to fulfill 
the assumed orthogonality between the line rj = const and the mirrors' worldlines, the function g(x,t) must satisfy 



fl (a. = 0,t)=0; g(x = L x (t),t) = 0; 
d x g{x = 0,t) = ; d x g(x = L x (t),t) = -L. 

There are many solutions to the above conditions, that can be written in the form 

g(x,t) = L x {t)L x (t)v(x/L x (t)), 



(16) 



(17) 



where i>(0) = u(l) = 0, v'(0) = and v'(l) = —1 (the prime denotes derivation with respect to x). For example, a 
possible solution is v(z) = \{z 2 ~ z A ). Physical quantities like the number of created particles or the energy density 
inside the cavity should not depend on the choice of g(x,t). We will keep a general function as a benchmark for our 
calculations. 

Finally, the mode functions in Eq.(|io|) can be expanded in terms of the instantaneous basis Eq.(^) as follows 



n (x, t>0)=J2 #>(,) J^r cos ( ) J± cos ( & ) J± cos ( ^ 



, Kv) \Kv) 

E^k n) (*)+^ n) (*)^.*)]^^ 

E^k n) W + QL n) W.9(^O]0k(x,L ;c W). 




k z ir 



kyir 



y l \ — cos 



(18) 



The initial conditions are given by 



/2uj t 



(0) 



(19) 



provided that L x {t) and L x (t) are continuous at t = 0, and that the initial acceleration L x (0) — 0(e 2 ). In this way 
we ensure that each field mode and its time derivative are also continuous functions at t = 0. 
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B. Dynamical equations 



We now study the trajectory L x (t) = Lq(1 +esin(Qt)). [] The equations for the modes Q^(t) can be obtained from 
Eg. (|i~8|) , since □w n (x, t > 0) = 0. We first apply the Dalambertian, and then multiply by </> k and integrate over the 
cavity. To order 0(e 2 ), the equations read 

Ql n) + Qi n) = -2A(t) £ g ik Q\ n) - X(t) £ g-^Qf 

j j 

-2X(t)Ll(t) £ r jk Qj n) - £ Qj n) (rjkA(t)L»(t) - A(^ jk ) - A(i)L|(t) £ r^Q^, (20) 
j j j 

where 

/ / i. \ 2 / ;„ \ 2 / /„ \ 2 

(21) 
(22) 
(23) 

dz [(»" - wfu)#j0 k + 2u>j<£ k ]. (24) 

Here, Uj is the frequency of the mode evaluated at e = 0. As before, the prime denotes derivation with respect to x. 
The coefficients gj k are defined by 




gjk = ix(i) / dx / dy / dz — 4- k = <^ v ' Wll fc * J » x -r jx, ( 2 5) 

Jo Jo Jo [_ - ofc y j y Ofc^ it k x = j x . 

Had we considered Dirichlet boundary conditions on the static walls (y — 0, L a ; z = 0, L z ) we would had obtained 
the same dynamical equations for the modes Q k (i); i.e., the form of the equations only depends on the boundary 
conditions imposed along the x-direction. This is because the coefficients rj k , rjij k and gj k do not depend on the 
particular form of the k in the plane y — z, as long as they are properly normalized in this plane. However, when 
Dirichlet boundary conditions on x — and x — L x (t) are considered, the equations for the modes are different (see 
Ref. [|l3| and Eqs.(^) and ([49]) below. Note that the coefficients c/j k for Neumann boundary conditions are different 
from those for Dirichlet boundary conditions. 

When the mirror returns to its initial position for t > ig na i the rhs of Eq. (^(i| ) vanishes and the solution reads 

Q k n) (i > tfi„ai) = 4 n) e^ kt + B^e-^, (26) 

with and being some constant coefficients to be determined by the continuity conditions at t = tgnai- The 
number of particles in the mode k is given by 

(AA k )=^2 Wk |4 n) | 2 . (27) 



1 Strictly speaking, we should add to L x (t) some decaying function in order to meet the continuity conditions at t = 0. Since 
we will be interested in a resonant behavior of the field, this additional function will not contribute, being irrelevant for what 
follows. For a more detailed discussion of this point see ttfX. 
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C. Multiple scale analysis 



In order to find a solution to Eq.(pO[) we use the multiple scale analysis technique Jig], which we have already 

£ Q k n)(1) . 

k , , „ -r^ n) (r)e _lWk *. The functions 



applied in our previous paper Jl3| . We first introduce a second timescale r — et and expand Q k n ^ = 
Replacing this into Eq.(pO|) we obtain, as zeroth order solution, — (r)e 4 " k * 



A^\t) and B^\t) will be obtained from imposing that no secularities appear in the equation for ■ This reads 



l(n)(l) 



-2L 2 n 2 sm{Qt)r kk ujlQ^ m + L 2 VL cos(fti)r kk u;£Q k n)(0) 



-L 2 ncos(nt)Q^ m [n 2 r kk + - 1 r ]kk \ 



+n 2 sin(nt) 5jkQj n)(0) - 2^ cos(nt) fjkQj 



i(")(0) 



j^k 



-2L 2 n 2 sin(nf) ^ r jk w|Qj 

j#k 



(n)(0) 



j#k 

Lgfi COs(ftt) rjfcWj $j 
j#k 



(n)(0) 



-L^cos(m)^Q 



(n)(0) 



j#k 



fi2? > + T2 r /Jk 



(28) 



where we have used that, to zeroth order, = — <-^ k Q k n ^ ^ ■ 

The equations for (r) and B^ (r) are obtained imposing the condition that any term in the right-hand side of 



Eq. (p8[) with a time dependency of the form e ±IWk * must vanish. We get 



dA 



(n) 



dr 



1 

2wk 

-E 

j#k 

E 

j#k 

E 

j^k 



£8 



2Ji 



9., 



B^ ] S{2uj k - Q) 



s(—uj k — ujj + n) 



2w k J 



(n) 



^yjy h zjz J 



5(w k -cjj + fi)-^A[ n) 



(29) 



and an analogous equation for n \ obtained by the interchange A^ «-> -B k "' ) . Note that Eq. (|29|) is independent of 
g(x,t). This non trivial check of our calculations follows from two identities we used to derive Eq. (f29[) , namely 



w k r jk 



2u k 



T^Jk 



L 2 cj 2 B k n) S(2io k - fi) 



i.(t) 



dy I dz(4> k v ) 



M*) 



dx 



dy 



L- 



dz\v' fafa + «(0k^j - k ^j)]' 



w k 5 k n) <5(2u; k - fl), 



(30) 
(31) 



where we have used the boundary conditions v'(0) = and v'(l) = —1. 



D. Examples 



Let us consider the "parametric resonant case" , in which the external frequency f2 is twice the frequency of an 
unperturbed mode k, Q — 2uj k . A second mode j will be coupled to the mode k iff \uj k ± Wj| = ft. We first assume 
this is not the case. Therefore the evolution equations become 



G 





1 


rfc2 2 




4 n) 


dr 


2wk 


I Ll 


dig* 


1 






A {n) 


dr 


2w k 


i Ll 



(32) 



It is easy to check from these equations that 

k 1 ° 

where ui 2 = u> 2 



(n) and R (n) 



k grow exponentially as e AjvT , with a rate Ajy = ^"(^k+^p)' 
It is interesting to compare this rate with that for Dirichlet conditions, which is given by 



Ad 



> 1. 



(33) 



For a given mode, the rate for Neumann boundary conditions is always bigger than the rate for Dirichlet conditions. 
Let us now assume the existence of one mode, say j, that satisfies = 3a->k and j y = k y , j z — k z . We obtain for 



4" } and 



dA^ 


1 




dr 




[ Ll 


dB™ 


1 




dr 


2w k 


I Ll 



-2u)i 



Bi 



A 



(n) 



1 

2w k 
1 

2uj k 



(-1) 
(-1) 



T 2 

n 2 w 2 
T 2 



Qui 



-6w£ 



A 



(n) 



B 



(n) 



(34) 



We also assume that the spectrum is such that the mode j is only coupled to the mode k. The equations for Aj and 



b\™^ are therefore 



dA\ 



(n) 



dr 



dB) 



(n) 



dr 



1 

6wk 
1 

6w k 



(_ 1 )*.+3« 



We write the system of equations in matrix form 

dv(r) 



= Mv(t), 



where 



v(r) 



Af\r) 



M 



1 

2wk 



fc^ 2 

Ll 


f 2^ 


4 n) 


k^ 2 
Ll 


f 2t4 


4-) 



(35) 



(36) 



a 6 

a 

c 

c 



(37) 



where a 



2u£ 



6^ 



, and c 



(-1) 



2^ 



system can be easily obtained after diagonalizing the matrix M.. The eigenvalues are given by 

A = -L(± a ± 
4w k 



•46c). 



The solution to this 



(38) 



We note that the exponential growth rate in the uncoupled case is given by A at = a/2wk- When two modes are 
coupled, the rate is given by the real part of the biggest eigenvalue in Eq. (|38|) . When a 2 + Abe < 0, the rate is half 
the one expected for the resonant mode when the coupling is neglected, as for Dirichlet boundary conditions jl3) . It 
is easy to show that this is the case if (— l) fca =+^ = +1. However, in the opposite case, be > 0, and the growth rate 
for coupled modes is bigger than \n = a/2oj-^ . 



A relevant case where two modes are coupled is the cubic cavity L x 
cavity frequency, 



L,, 



L z = L, We fix SI as twice the lowest 
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n = 2 W(l!l)1) = (39) 

The fundamental mode k = (1, 1, 1) is coupled to j = (5, 1, 1) because W(5 j. i) = i Only these two modes are 

coupled, since there does not exist in the spectrum any mode s satisfying lu s — 5u)ni t i\. For this particular case, the 
four eigenvalues are 

A = — ^^(±5 ± 6.35i). (40) 

Had we neglected the intermode coupling, we would have concluded that the growth rate in the fundamental mode 
would be A = 2.5n/^/3L. The growth rate in the coupled case is one half of this. 

One striking new feature is the possibility to enhance the exponential growth rate by means of mode coupling, 
provided that the two coupled modes satisfy (— \) k * + iy = —1. As an example let us consider a cavity with dimensions 
L y = L z — 4L X . We set the external frequency to be 



2tt 



n = 2w (0) i il) = -^-. (41) 



If this is the case, then the mode k = (0, 1, 1) is coupled to j = (1, 1, 1). The four eigenvalues are 



A=^ z -(±l±v/3i/3). (42) 

This means that the exponential growth is at the rate 0.7Air/ L x , which is more than twice the value we would had 
predicted had we neglected the coupling (n/y/8L x ). 



IV. THE ELECTROMAGNETIC FIELD 



A. Transverse electric modes 



For the TE case, the expansion of the vector potential for an arbitrary moment of time, in terms of creation and 
annihilation operators, can be written as 

A( TE )(x,t)=^a>( TE) (x,t) + H.c. (43) 

11 

For t < the cavity is static, and each field mode is given by 



u( 1 TE) (x,t<0)= ' 



L X LyL Z 



0,a n sin f^r^J cos yj^vj sin (""7~ Z J >/3 n sin (jj^-xj sin (j^v) cos (y^ z )) e ^ 

where n x , n y , and n z are integers such that n x > 1, n y ,n z > 0, and n y ,n z cannot be simultaneously zero. The con- 
stants a n and /3 n are components of the polarization vector for the electromagnetic field, satisfying the normalization 
condition ct£ + 0^ = 1, and the Coulomb gauge condition, a n n y /L y + I3 n n z /L z = 0. 

When t > 0, we expand the mode functions in Eq.(^3|) with respect to an instantaneous basis 



uL TE) (x,t > 0) = EQgr^^y <*» *1™W)> ( 45 ) 

where <&[ TE ' ) is 
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The functions ^j^™" 1 form a complete set satisfying 

LV dyf* dz *£g • = 4„ j, ^ ■ (47) 

From the above Eq.(f45|) it is easy to obtain the dynamical equations for the modes QkTB- We S e ^ 

Qk/TE + ^(t) Qk°TE = 2A(i) 2 SmQ$E + A(f) £ <7kjQ$ E (i) (48) 

j j 

where 

flkj = -ffjk = n ^ y ' B * J * ., , _ (49) 

As expected, these equations are exactly those corresponding to a scalar field satisfying Dirichlct boundary conditions 
on the surfaces x = 0,L x (t) |Q. Therefore, the number of created photons in the TE mode equals the number of 
created Dirichlet scalar particles. 

As an example, let us consider the parametric resonant case f2 = 2wk for a cubic cavity. For uncoupled k modes 
(such as either of the two fundamental TE modes, k = (1, 1, 0) and k = (1, 0, 1)), the number of TE photons grows 
exponentially as 

(A/k,TE> =sinh 2 (A B et), (50) 

where Ad is the growth rate for Dirichlet scalar particles, introduced in Section III.D. For the above mentioned 
fundamental modes, Xd = n/2^/2L. The first coupled TE mode is k = (1, 1, 1), which only couples to the TE mode 
j = (5, 1, 1). At large times et/L^> 1 the number of TE photons in those modes grows as (A/k te) ~ (A/j te) ~ e°- 9e */ L 



B. Transverse magnetic modes 

The expansion in terms of creation and annihilation operators is again of the form Eq. ( ]43| ) , but now for t < each 
field mode is given by 

\/2uj n ^ L x L y L z 

^ ,a n cos x) sm yj cos zj ,/? n cos xj cos [j- yj sin zj je . (51) 

Here n x , n y , and n z are non-negative integers, and n y and n z cannot be simultaneously zero. 

On the other hand, when t > we introduce an instantaneous basis similar to that of the scalar field in Section III. 
We write 




where the functions *&L k are similar to their TE counterparts (they can be obtained by interchanging cos and sin in 
the RHS of Eq.([46|). Since all TE modes have n x = 0, the first mode of the cavity that can be excited by the external 
frequency is a TM mode. In particular, for a cavity such that L x <C L y , L z only TM modes can be excited. 

From the above equation, it is now clear that the dynamical evolution of the TM modes is that of a scalar field 
satisfying generalized Neumann boundary conditions. As a consequence, the number of created photons in the TM 
mode equals the number of created Neumann scalar particles. Again, let us consider the parametric resonant case 
fl = 2w k for a cubic cavity. For uncoupled k modes (such as either of the two fundamental TM modes, k = (0, 1, 0) 
and k = (0, 0, 1)), the number of TM photons grows exponentially as 

(A/k,TM> =sinh 2 (A iV et), (53) 
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where Xpf is the growth rate for Neumann scalar particles, also introduced in Section III.D. For the fundamental 
modes A^v = ir/L. The first coupled TM mode is k = (0, 1, 1), which only couples to the TM mode j = (4, 1, 1). For 
large times (et/L ^> 1) the number of particles in these modes grows as (A/Ic.tm) ~ (A^,tm) ~ e 4 4et / L . The next 
coupled TM mode is the same as the TE mode, namely k = (1, 1, 1), coupled to j = (5, 1, 1). The exponential growth 
is (Mc.tm) ~ (A/j,TM) ~ e 4 - 5et / L , the growth rate for these modes being greater than that for the TE case. 

V. DISCUSSION 

In this paper we have computed the resonant photon creation inside a three dimensional oscillating cavity taking the 
vector nature of the electromagnetic field into account. Previous works studied the case of a scalar field with Dirichlet 
boundary conditions. As the electromagnetic field involves both Neumann and Dirichlet boundary conditions, we first 
analyzed a massless scalar field satisfying (generalized) Neumann boundary conditions. We have shown that in this 
case it is also possible to expand the field modes in terms of an instantaneous basis, the difference with the Dirichlet 
case being that the expansion is not unique - it depends on an arbitrary function g(x,t) satisfying the boundary 
conditions Eq.(|l6|). However, physical quantities like the number of created particles or the energy density inside the 
cavity are independent of the choice of such a function. After treating the Neumann scalar case we considered the 
full electromagnetic problem and showed that the TE modes of the elecromagnetic field are essentially described by 
a Dirichlet scalar field, while the TM modes correspond to a Neumann scalar field. 

We have studied in detail the resonant situation £1 = 2u>k for two cases: an uncoupled resonant mode and two 
coupled resonant modes. In both cases, the exponential growth of created photons is greater for TM modes. For the 
uncoupled case, we have found that 

a™ = ^4 

Ate ~ w p 

For a cavity with L x ~ L v ~ L z , uj^ ~ |w£ so Atm — 5 Ate- We can estimate the number of created TE and TM 
photons given by Eqs.( p0| , p3| ) using typical values for the maximal dimensionless displacement e that may be obtained 
in conceivable future experiments. For 3D cubic cavities of linear dimensions of the order of 1 — 10 cm, the lowest 
resonant frequency is of the order of GHz. It may turn out to be very difficult, if not impossible, to make the cavity 
oscillate as a whole at such a high frequency. To overcome this difficulty a different experimental scenario was proposed 
in p2| , consisting of strong acoustic waves excited on the surface of the cavity wall. Typical materials cannot bear 
relative amplitude deformations in excess of <5 max = 10~ 2 . This sets a limit to the maximum velocity of the boundary, 
Wmax = imm" s ~ 50m/s, (v s is the speed of sound in the material), and consequently to the maximal dimensionless 
displacement e max = v max /f2L. For example, for a cavity with L = 10cm whose lowest mode (i.e., either of the two 
TM modes k = (0, 1, 0) or k = (0, 0, 1)) is being excited (ft = 2nc/L = 18GHz), we get e max « 10~ 8 . Even for a value 
of e 10 times smaller than this, one gets an exponentially large number of created photons (.Mc.tm) = sinh 2 (10i/s) 
which, after 1 second, gives a total of approximately 10 8 photons created in that mode. We can also compare the 
number of photons produced for an uncoupled mode k, common to both kind of polarizations TE and TM. For 
example, for the mode k = (1, 1,0) one gets (Mc.te) ~ sinh 2 (3i/s) and (Mc.tm) ~ sinh 2 (10£/s), which after 1 second 
produces a total of 10 2 TE photons and 10 8 TM photons. For the case of two coupled modes we have found that, 
for Neumann boundary conditions, the coupling can enhance the exponential growth. This contrasts with the case of 
Dirichlet boundary conditions, in which the coupling always suppresses the exponential growth. These facts may be 
relevant for an eventual experimental verification of the dynamical Casimir effect. 

All the above considerations assume ideal conditions, such as perfectly conducting plates, exact parametric resonant 
condition il = 2^, arbitrary large Q factor for the cavity (no leakage of photons), no thermal noise, etc. Some of 
these conditions were relaxed in our previous paper fis}] , where we analyzed, for Dirichlet boundary conditions, the 
enhancement of photon creation due to finite temperature effects, slightly off-resonance situations, the case of three 
coupled modes, etc. The generalization of these findings to the electromagnetic case is not too complicated, and we 
expect similar conclusions. Given our results for the dynamical behavior of TE and TM modes, it is also possible to 
study the full electromagnetic problem in three dimensional leaky cavities along the lines of [flSj] . 
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